Canonical quantization of general relativity in discrete space-times 
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It has long been recognized that lattice gauge theory formulations, when applied to general rela- 
tivity, conflict with the invariance of the theory under diffeomorphisms. Additionally, the traditional 
lattice field theory approach consists in fixing the gauge in a Euclidean action, which does not ap- 
pear appropriate for general relativity. We analyze discrete lattice general relativity and develop a 
canonical formalism that allows to treat constrained theories in Lorentzian signature space-times. 
The presence of the lattice introduces a "dynamical gauge" fixing that makes the quantization of 
the theories conceptually clear, albeit computationally involved. Among other issues the problem 
of a consistent algebra of constraints is automatically solved in our approach. The approach works 
successfully in other field theories as well, including topological theories like BF theory. We discuss 
a simple cosmological application that exhibits the quantum elimination of the singularity at the 
big bang. 



Lattice approaches have been successful in Yang-Mills 
theory not only in making the theory finite, but also in 
making it practically computable. There has been for 
a long time a hope that similar techniques could help 
in defining a quantization for general relativity (see jj[ 
for a review). More recently, discretizations on lattices 
have been used to regularize the Hamiltonian constraint 
of general relativity by Thiemann (2) and others and also 
to define the path integral in the "spin foam" approaches 
(see H for a recent review). A problem that has been 
present since the outset is that the introduction of a 
lattice breaks the symmetry of general relativity under 
space-time diffeomorphisms. This has had several prac- 
tical implications. For instance, discretizations of the 
constraints of general relativity fail to close an algebra 
like the one present in the continuum Q. Moreover, it 
is well known that if one discretizes the Einstein equa- 
tions, the evolution equations fail to preserve the con- 
straints. Therefore, strictly speaking, the discrete theory 
is inconsistent. This has been noted and discussed in the 
numerical relativity literature 

In this letter we would like to point out that what has 
been lacking is a systematic canonical treatment of the 
discrete theory, both at a classical and quantum mechan- 
ical level. When one treats the discrete theory carefully, 
one notices that it has an involved canonical structure. 
Dealing with it is non-trivial, since one needs a general- 
ization of the Dirac procedure to systems in which time 
is discrete (attempting a purely spatial discretization ap- 
pears highly unnatural in the case of general relativity, 
where the constraints mix space and time symmetries). 
We present such a generalization of the Dirac proce- 
dure and show that the resulting theory is consistent, 
in the sense that constraints do not conflict with evolu- 
tion equations and can be quantized. In our construction 
the lattice provides a natural "gauge fixing" for the the- 
ory, which then has a "true Hamiltonian" evolution. This 
gauge fixing may have practical use in classical numerical 



relativity as well. 

Let us illustrate our procedure with a mechanical sys- 
tem to begin with. We take time to be discrete. The 
action therefore is written as S = XL,=i -^(ini 1n+i) 
where we have partitioned the time interval and dis- 
cretized the derivatives of the Lagrangian. The Lagrange 
equations of motion are, dS/dq n = dL(q n _x,q n )/dq n + 
dL(q n ,q n+ i)/dq n = 0. In order to introduce canon- 
ically conjugate momenta one notices that differenti- 
ating with respect to the time derivative is equiva- 
lent to differentiating with respect to q n +i, or, alter- 
natively, p n = dL(q n -i, q n )/dq n . The evolution of 
the canonically conjugate pair (q n ,Pn) can be induced 
through a type 1 canonical transformation generated by, 
F\(q n ,q n+ i) = —L(q n ,q n+ i). The resulting equation of 
motion is p n +i = dL(q n , q n+ i)/dq n +i, and together with 
the definition of the momentum completely determines 
the evolution of the system and reproduce in the limit 
the equations of motion of the system j(| . 

To quantize the resulting system we need to find a uni- 
tary operator that implements the finite canonical trans- 
formation. Let us particularize to the example of a par- 
ticle in a potential. In such case the discrete Lagrangian 
is L(n, n + 1) = m{q n +i - q n ) 2 /{2e) ~ V{q n )e, where e 
is the time interval. The canonical transformation yields 
equations of motion that can be rewritten as determin- 
ing q n+ i,p n+ i in terms of q n ,p n . In order to quantize we 
need to find a unitary operator that implements quantum 
mechanically the equations of motion, p n +i — Up n U* , 
q n +i = Uq n W, This quantization ensures that q n ,p n 
are such that [q-mPn] — ih- The operator can be read- 
ily found. Starting with wave functions that are func- 
tion of the coordinates *£[q], the operator is given by, 

tf = cxp(^)exp(fe). 

We therefore see that one can set up a canonical theory 
for a system with discrete time and quantize the system. 
What we need to do in order to generalize this to gravity 
is to consider systems with constraints. The construction 
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proceeds as before, but the resulting discrete equations 
are radically different from the continuum ones: they 
are generically inconsistent unless one chooses particu- 
lar values of the Lagrange multipliers. These particular 
values also make the constraints automatically satisfied: 
the resulting theory is only described through evolution 
equations. As before, the quantization of the theory is 
achieved by finding a unitary transformation that imple- 
ments the finite canonical transformation. This is readily 
done for the case of the particle, at least for simple po- 
tentials. Notice that once we solved the constraints for 
the multipliers, they do not play any role in the quan- 
tization procedure. T. D. Lee f?J was the first to rec- 
ognize that the discretized mechanics of a particle in a 
potential could be made consistent through fixing the 
Lagrange multipliers (in his language, fixing the time in- 
terval). Also Friedmann and Jack noticed similar things 
in the context of Hamiltonian Regge calculus || . 

Solving the discrete constraints by choosing the La- 
grange multipliers is a procedure that is only available in 
the discrete theory. This procedure might be useful in 
numerical relativity. There it is well known that the dis- 
crete evolution equations (for instance, in the ADM for- 
mulation) fail to preserve the discrete constraints. Our 
proposal here consists in the following: take an evolu- 
tion step with arbitrary, unspecified, lapse and shifts. 
Then impose the discrete constraints by solving the for 
the lapse and shift. There are four quantities to be de- 
termined and four equations. The construction obviously 
is only guaranteed to work locally both spatially and in 
time. Numerical experimentation will be needed to de- 
cide if this is really useful in practice. 

We have tested the canonical procedure with Yang- 
Mills theories and BF theories in the lattice and one re- 
covers in the former case the traditional lattice theory. 
To our knowledge, our approach provides the first lattice 
treatment for BF theory ever presented. Details can be 
seen in ||. 

Let us now sketch the case of general relativity. We 
start by considering the Lagrangian written in terms of 
Ashtekar variables, (see for instance |j) and the book by 
Ashtekar (To) page 47). To simplify treatment here we 
will consider the Euclidean case, but there is no obstruc- 
tion to treating the Lorentzian case, one simply has to 
replace the action we are considering here with the one 
introduced by Hoist JTTJ which leads to the real variables 
of Barbero jl2). The Lagrangian is, 

L = J E ai F l a0 + e abc [E u E ci e ijh N + N b E ck ] F$ e e ade , 

where E al and A J a are the usual Ashtekar variables Jlj| 
and N and N a are the lapse and shift. This Lagrangian 
can be discretized as follows, 

L(n,n+ 1) = £Tr [E^ v h a n>v V n , v+ei (h a n+h J(V n , v y 



+ K l,n,vh 2 n,v h3 n,v+e 2 ( h2 n,v+e 3 )Hh 3 n ,y + CJC. 

+ a a , n , v (h a ntV 0^) f - l) + p ntV (V n , v Vl v - 1)] (1) 

where K a ^^ v = e abc [E^ v E^ v N n , v + N^ V E^ V ] and re- 
peated indices a, b, c are summed from 1 to 3 and "eye" 
means cyclic permutations of 1,2,3. In this expression 
hj l+1 represents an holonomy along the a direction at in- 
stant n+l,V n represents the "vertical" (time-like) holon- 
omy. We will assume that the holonomies are matrices 
of the form h = h^T 1 , V = V'T 1 where T° = I/V2 and 
T a = -ia a /V2 where a a are the Pauli matrices. The 
indices n, v represent a label for "time" n and a spatial 
label for the vertices of the lattice v. The elementary 
unit vectors along the spatial directions are labeled as 
e a , so for instance n-\-e\ labels the nearest neighbor to n 
along the e\ direction. The quantities E% v are elements 
of the algebra of su{2) and a and (3 are Lagrange multi- 
pliers, the last two terms of the Lagrangian enforcing the 
condition that the holonomies arc elements of SU (2) 

We need to work out the equations of 
motion of the action. The variables are 

^,«. ^,». j3^„, JV ni „, JV^„. For each vari- 

able we compute the canonical momenta at instants n 
and n + 1 by dL/dq n and dL/dq n +i as in the examples 
we discussed. We will not list them explicitly for reasons 
of space but we will discuss their implications. The 
equations for a lead to a constraint that implies the 
spatial holonomy h is a matrix of SU(2). The equations 
for (3 lead to a similar conclusion for the holonomy V. 
The equations for E lead to an evolution equation for the 
holonomy. The equations for h lead to a constraint relat- 
ing E with the canonical momentum of the holonomy at 
instant n + 1 and can be viewed as an evolution equation 
for E. The equations for N, N a give the diffeomorphism 
and Hamiltonian constraints of general relativity, which 
are solved to determine the multipliers N, N a as we 
discussed above. The equations for V yield Gauss' law. 
To complete the canonical theory one needs to check 
that all constraints are preserved in time. This can be 
done and it leads to a set of second class constraint 
very similar to the ones we have analyzed in detail in 
reference Q for BF theory. One needs to introduce the 
Dirac brackets, which are identical to those of BF theory 
and impose the constraints strongly. This concludes the 
canonical formulation. 

Quantization would now proceed by finding a unitary 
evolution operator that implements the previous evolu- 
tion equations as operator equations quantum mechan- 
ically. Notice that there is no issue of algebra of con- 
straints, since the latter have been solved for the multi- 
pliers, the lapse and shift. The procedure is computa- 
tionally intensive in a generic situation, requiring the si- 
multaneous solution of four non-linear coupled algebraic 
equations at each lattice point. For particular situations, 
the procedure might be completed in a straightforward 
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fashion. The study of midi-superspace situations is there- 
fore the next natural step in this program of quantiza- 
tion of gravity. A direct application would be to study 
the dispersion of waves to attempt to make contact with 



gamma-ray-burst phenomenology 1 14 

For concreteness, let us exhibit the construction in a 
simplified situation, that of a cosmological model. Wc 
consider a Friedmann model with flat spatial sections 
coupled to a scalar field which we take to have a large 
mass so we can neglect its kinetic term in the action. 
This is just done in order to yield the simplest possible 
cosmological model with at least one continuum degree 
of freedom. Written in terms of Ashtekar's variables, the 
only non-trivial quantities are one component of the triad 
field E(t) and one component of the connection A{t) and 
for the scalar field its value <f>(t) and its canonically con- 
jugate momentum II(<). The continuum Lagrangian is 
given by L = EA + U<j> - NE 2 (-A 2 + AE + m 2 <j) 2 E) 
where N is the lapse (with density weight — 1 , as is usual 
in the Ashtekar formulation, in terms of the usual lapse a 
it is given by A = E~ 3 ^ 2 a) and A is a (positive) cosmo- 
logical constant. The exact solution for this model (in the 
time slicing in which a = 1) has a constant scalar field 
and the geometry is given by (anti)DeSitter space, E — 
exp(2A/A + m 2 cj) 2 t)/(A + m 2 (f> 2 ), A = exp(yA~ 
we discretize the action we get, 



<j) 2 i). If 



L n = E n (An+1 - A n ) + 7T„ (<j) n +l - 4>n) 

-N n E 2 n (-A 2 n + AE n + m 2 <\> n 2 E n ) 



(2) 



and the generating function of the canonical transforma- 
tion that implements the equations of motion is given as 
usual by F\ = —L. Working out the equations of motion 
and the definition of the canonically conjugate momenta 
from the generating function and the Lagrangian, one 
ends up with a recurrence relation for the variable A n . 



A 



A 2 n + 2 (A n+1 - An) An = 



(3) 



The lapse is completely determined by solving the 
Hamiltonian constraint. This fixes dynamically the 
gauge: the evolution step in time is determined by the 
initial conditions. For instance if we substitute the ex- 
pression for the lapse in the equations of motion one 
gets, E n = P A n +i = A n 2 /(A + m 2 <f> n 2 ), A n+1 = A n + 
(Al - AP A n - m 2 (t> n 2 P A n ) /2A n where P A is conjugate 
momentum to A. We see that for instance the "step" be- 
tween A n and A n+ i is given by A 2 n — AP A n — m 2 (j> n 2 . 
Choosing the initial values for A and P A and the scalar 
field determines the step. 

This simple model exhibits several attractive features. 
For instance, for large values of the "time label" n, 
A -> n 2 / 3 , E -> n 4 / 3 and the lapse N -> n~ 3 . If 
we identify t — en where e is a constant, we then get 
E — > £ 4 / 3 for large times. This approximates the con- 
tinuum solution if one chooses a lapse in the continuum 
a r~j t . We can recover any parameterization we desire 



in the continuum by redefining t as a function of en. We 
therefore see in which sense the lattice "fixes the gauge" . 
Indeed, the lattice treatment operates just like a gauge 
fixed treatment, but one still has available the full repa- 
rameterization invariance of the theory when identifying 
the lattice behavior with the continuum one. Although 
not a traditional gauge fixing, this procedure may share 
with the latter some problems. One breaks the gauge 
symmetries of the theory, but it is yet to be seen, espe- 
cially in complex situations if one does not face a Gribov 
problem or that singularities do not appear that prevent 
from covering portions of the gauge orbits, etc. In par- 
ticular, by choosing different parameterizations we can 
choose to start the region that approximates the contin- 
uum behavior as early or as late as we desire. One can 
consider the recursion relation ([}]) and run it backwards 
in time. If one does that, one finds that the universe 
achieves a minimum radius and the recursion continues 
with the universe re-expanding again. This is a dramat- 
ically different behavior than the one in the continuum 
theory. In the latter, if one gives any initial condition and 
runs back in time one will always meet an initial singu- 
larity (in this simple model, just a coordinate singularity, 
in other models it is a genuine singularity [pp[ ). In the 
discrete theory one needs to fine tune the initial data for 
this to happen, otherwise the singularity is avoided: it is 
a point that does not fall on the grid. Quantum mechan- 
ically, the implications of this fact are more significant. 
Since the singularity is only achieved in a set of measure 
zero of the classical theory, quantum mechanically one 
has zero probability of hitting the singularity. Notice 
that this mechanism for avoiding the singularity quan- 
tum mechanically appears distinct, although it predicts 
somewhat similar results, to the one recently presented 
by Bojowald [flEf in the context of loop quantization. Sin- 
gularity avoidance has also been discussed in other con- 
texts JiJI with different mechanisms. We would like to 
stress that the particular way the singularity is avoided 
is model dependent and in some models the singularity 
may not be avoided |p0| 

Quantization proceeds by considering wave functions, 
for instance, of the momentum (P — P ) and the scalar 
field, Py{P,cj>) = P^(P,4>), A$(P,4>) = <N>(P,4>)/dP. 
One needs to implement the reality conditions A* = —A, 
P* = P. These can be implemented via the inner prod- 
uct < *|$ >= J^ o d0/^ o dF**(F,(/))$(P,0). Since 
under this inner product P has a continuous spectrum, 
the probability of having a singularity (which corre- 
sponds to P = 0) vanishes. This can also be seen from 
the fact that the expectation value of the volume is al- 
ways greater than zero. This implements in detail the 
behavior we had anticipated in the heuristic discussion 
above. 

Going back to the generic discussion of our approach, 
at this point one might wonder if by taking too seriously 
the discrete theory we are not neglecting the continuum 
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limit. Two comments can be made. First of all, there 
appears to be some consensus among researchers in 
quantum gravity that the ultimate theory may have some 
fundamental discreteness at the quantum level. This is in 
part based on the fact that the continuum limit of Yang- 
Mills theories relies heavily on the renormalizability of 
the theory, which appears not to be present in the gravi- 
tational case. This is therefore motivation to take a look 
at the discrete theory seriously. Moreover, having a con- 
sistent discrete theory, which includes among its solutions 
some that correspond to solutions of the continuum the- 
ory, opens the possibility of defining an averaging proce- 
dure that could make results discretization-independent. 
If one adds up solutions for various discretizations, the so- 
lutions that have a continuum counterpart will appear in 
all terms and will add up. The solutions without a con- 
tinuum counterpart will therefore be suppressed, since 
they will be different for different discretizations. No- 
tice that this allows a method of implementing the "sum 
over all discretizations" that is advocated in the context 
of spin foam models []l~7| . Following the procedure out- 
lined above for gravity, one can construct transition am- 
plitudes for the discrete theory. Then, performing the av- 
eraging we just discussed for a finite, yet large, number of 
discretizations will produce a result for the transition am- 
plitude that is approximately discretization-independent 
(this is only a heuristic argument that could face prob- 
lems, for instance, discretizations could be unstable and 
modes that do not correspond to the continuum limit 
could end up dominating the sum). If the construction 
were to work one would not need renormalizability in or- 
der to define a continuum theory in this way. It should be 
noted that the general state of the art of these ideas is, at 
the moment, only speculative, and there is disagreement 
among various researchers. In fact, the more conven- 
tional viewpoint for the continuum limit of lattice gauge 
theories, applied to quantum gravity does appear to work 
in sucessfully in I + 1 dimensions, and progress is being 
made in higher dimensionalities JT^]. Even within this 
more traditional viewpoint, the "consistent" construction 
we present in this paper could play a role at the time of 
selecting the correct measure of integration, since it au- 
tomatically takes into account the constraints, which is 
crucial at the time of properly defining the path integral. 

Another attractive element of this treatment is that 
the separation of the lattice points is dynamically deter- 
mined by the initial conditions. This implies that quan- 
tum mechanically, where generically one has a weighted 
superposition of all possible initial conditions for the sys- 
tem, there is an automatic "averaging" of all possible dis- 
cretizations per each given quantum state. It should be 
noted that this averaging feature only occurs for totally 
constrained systems so it is a case that diffcomorphism 
invariance actually helps in getting rid of the details of 
the particular discretization. 

Finally, the lattice treatment introduced for gravity, 



since it operates as a gauge fixing, provides a solution for 
the "problem of time" . However, as it happens with usual 
gauge fixings one should expect that our construction will 
only work locally. 

Summarizing, we have presented a formulation of 
general relativity on the lattice that is consistent and 
well defined classically and can be readily quantized. 
The formulation can operate in both the Euclidean and 
Lorentzian signatures and yields naturally a proposal for 
the continuum limit of the theory. We exhibited its be- 
havior in detail in a particular cosmological situation 
where one sees the singularity at the big bang disappears 
due to quantum effects (although this is not a generic 
feature, see Q). The immediate future course of ac- 
tion is to explore the formulation in more realistic, midi- 
superspace models, where detailed predictions of the field 
theory aspects of quantum gravity could be worked out, 
some of which may have experimental implications. 
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